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Abstract 

A quantum statistical theory of spin-dependent tunneling through asymmetric magnetic double 
barrier junctions is presented which describes both ballistic and diffuse tunneling by a single an- 
alytical expression. It is evidenced that the key parameter for the transition between these two 
tunneling regimes is the electron scattering. For these junctions a strong asymmetric behaviour in 
the I-V characteristics and the tunnel magnetoresistance (TMR) is predicted which can be con- 
trolled by an applied magnetic field. This phenomenon relates to the quantum well states in the 
middle metallic layer. The corresponding resonances in the current and the TMR are drastically 
phase shifted under positive and negative voltage. 



PACS numbers: 72.10-d, 72.25-b, 73.40Gk 



After the discovery of a large room temperature magnetoresistance effect in magnetic 
tunnel junctions (MTJ) |]J many potential applications have emerged based on the spin 
polarized transport through a thin insulating barrier sandwiched between two ferromagnetic 
metals|I], ^, ^ ffl. The current state of the art production of tunnel junction elements with 
a resistance variation of up to 40% when magnetically switched []5], |5| makes them very 
promising in particular for the application as tunneling magnetic random access memories 
(MRAM). Current MRAM designs || [7| which incorporate arrays of magnetic tunnel junc- 
tions add an additional semiconductor switch in series (either a CMOS transistor or a p-n 
junction) with the MTJ memory cell to suppress (or block) parasitic signal paths within the 
array of lines. However, such a concept is hampered by the difficulty of mixing semicon- 
ductor and metal technology. In this letter, we propose a novel theoretical concept which 
eliminates the introduction of additional semiconductor components by using an asymmetric 
double barrier structure M1/O2 a/M 3 6/O4 c/M 5 (Mj and Oj are the magnetic metal and 
oxide layer respectively with corresponding thickness a, b and c) that in itself acts as the 
blocking device. It is shown that the I-V characteristics and the TMR of such asymmetric 
double barrier structures have a strong "diode" -like behaviour under positive and negative 
applied voltage. Moreover, the asymmetric properties can be varied by a magnetic field 
leading to the concept of a "magnetically controlled" diode. 

The potential profile of the system under applied voltage V ext is shown in Fig. 1(a), where 
Ui and Vj (i=2,4) are respectively the potential of the barrier region and the linear voltage 
drop therein, and V° (i=l,3,5) is the spin-dependent potential of the z-th metal. The outer 
metallic layers are assumed to be semi-infinite. The symmetric structures with a = c and 
U 2 = U4 were considered in references |j) and |J but the resonance tunneling was described 
without taking into account the electron's scattering inside the middle metallic layer. In this 
case it was shown that the conductivity and the TMR exhibit resonance peaks as a function 
of the thickness of the middle metallic layer due to quantum well states. Here we will 
show that the asymmetry of the structure in combination with the presence of quantum well 
states leads to a large asymmetry in the current for forward (positive) and reverse (negative) 
applied voltage. The evaluation of the current through the double barrier junction is based 



on the determination of its transmission probability D pi. Here, D is derived using the 



Green functions technique in the mixed real space-momentum representation jnj |T3, 13 
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FIG. 1: (a) Potential energy diagram and (b) the calculated current-voltage curve for the asym- 
metric structure Cu/0 2 7 A/Cu 5.5 A/0 4 21 A/Cu with U 2 - Ep = U 4 - Ep = 3 eV. 



and it can be written in the form: 
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where G K = G K (z,z') and G* K = G* K (z,z') are retarded and advanced Green functions of 
the system, z is the coordinate perpendicular to the plane of the structure and k indicates 
the in-plane momentum of the electron. This expression is similar to the Kubo formula 
for the non-local conductivity [0 by replacing Fermi energy Ep with an arbitrary energy 
E over which an integration has to be done in the limits from Ep upto Ep + eV ex t- It is 
convenient to use eq. (|l|) since for the calculation of the Green function of a system including 
scattering it is possible to use well defined quantum statistical methods. In the presence of 
scattering in the metallic layers, the problem can not be solved exactly and the tunneling 
through the whole system is described by consecutive tunneling through each barrier. In 
this case the current densities in the first and second barrier are calculated separately and 
the condition of constant current throughout the structure has to be fulfilled by introducing 
either an effective electrical field Hl3| , [14| inside each barrier or by calculating so called vertex 
correction |1~3|, |l5|j. Here the first approach is used. The Green function is the solution to the 



Schrodinger equation in each layer, where in the barrier region the WKB approximation [16 
is used: 
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and 



where k a Ei = ^(2 mi /h 2 )(E - Vf), If and rrii are the mean free path of the electron with 
energy E and spin a and its effective mass in the z'-th layer, respectively, and voltage 

z < 

(zV 2 )/a < z < a 

Vi(z) = <V 2 a < z < a + b 

V2 + (z — a — b) V4/ c a + 6<^<a + 6 + c 
V e:rt 2; > a + b + c 

The solutions in each layer are matched by the boundary condition at the interfaces z — 0, 
a, a + 6, a + 6 + c (see Fig. 1(a)). From the Green function the non-local probability D 
is calculated using the formula ([!]). With the assumption that the voltage drop occurs 
only across the barrier regions (V2 and V4 in Fig. 1(a)) and that the conduction band 
edge is flat in the metallic regions, it is sufficient to consider only the non-local probability 
for those z and z' lying in the barrier regions 2, 4 0] (see Fig. 1(a)). This yields the 
four analytical expressions: D 22 , -D24, -D42 and _D 44 . It is found that these quantities are 
independent of z and z' and hence the condition that the divergence of the current density 
is zero is automatically satisfied inside the different regions. With these two point tunneling 
probabilities the current density throughout the first (j'2) and the second (j'4) barrier can be 



calculated separately and are written in the following form ref. [|10| 

#(4) = ^ / dE[f(E) - f(E + eV 2 )] J D^ m {E, k)kcIk 



J dE[f(E + eV 2 ) - f(E + eV ext )} J D° m) (E, k)k<Ik 



(2) 



where f(E) is the thermal occupation probability of a state with energy E. This is the gen- 
eralization of the linear response in metals where the current density at a point z is related 
to the electric field at a point z' through the two-point conductivity a(z,z'). Finally, from 
the requirement that the current density has to be constant throughout the whole barrier 
structure, j% = j%, the effective electric field inside each barrier and hence the voltage drops 
V2 and V4 are determined in a self consistent way for a given applied voltage V ext . The 
additional condition of the current conservation here results from introducing the effective 




FIG. 2: Dependence of (a) the forward and reverse current and (b) the corresponding asymmetry 
ratio on the layer thickness b of the middle metallic layer for the non-magnetic asymmetric double 
barrier structure (see text). 



electric field instead of calculating the vertex corrections J14| . Alternatively, if one calcu- 
lates the current taking into account the vertex correction, the matching of the boundary 
conditions would be enough to provide the current conservation. The resulting dependence 
of the current density on V ex t is given in Fig. 1(b) for the case of an asymmetric double 
barrier structure a ^ c or U 2 7^ U4, revealing a strong asymmetry in the I-V characteristic 
reminiscent of a diode. To understand the physical origin of this enhanced asymmetry, we 
first consider the case where Mi, M3 and M5 are non- magnetic metals, using for the Fermi 
wave vector the value of Cu k,E F i = 1.36 A -1 fpLTjl - Cu is chosen here because of its simple 
electronic structure similar to free electrons. Furthermore, we use the electron mean free 
path of 100 A for the scattering in the metallic layers and an electron effective mass of 
^2(4) =0.4 [0 inside the barriers. In Fig. 2(a) the current density is shown as a function 
of the thickness b of the middle layer M3 for a fixed applied voltage V ex t = 0.7 V for a 
double barrier structure of Cu/C>2 21 A/Cu b A/O4 21 A/Cu with different barrier heights 
(U 2 — Ep = 1 eV and U4 — E F = 3 eV). Both forward (positive V ex t) and reverse (negative 
V ex t) currents exhibit resonance peaks (oscillations) which are associated with the formation 
of quantum well states in the middle metal layer 3 |§, [El], [16], [0| [2(| [21], |2"2fl . The period 
of these oscillations is the same for both curves and proportional to ir/kE F (see, for example, 
ref. [[ill), but the positions (phases) of the resonant peaks are shifted with respect to each 
other. The phase is defined by the boundary conditions at the metal/oxide interfaces. In the 
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FIG. 3: Dependence of (a) the TMR ratio and (b) its corresponding asymmetry ratio on b calculated 
from the forward and reverse current through an asymmetric double barrier structure with magnetic 
outer layers (see text). 



case of asymmetric double barrier structures, the matching of the phases at the interfaces 
is sensitive to the direction of the current due to different D 2 2 and -D 44 in (g) which depend 
exponentially on the barrier parameters leading to asymmetric voltage drops V2 and V4. 
The difference in these voltage drops will bring the quantum well states in M3 to line up 
differently with respect to the energy E of the electron under positive and negative applied 
voltage. The resulting phase shift of the current density leads, thus, to a current asymmetry 
ratio (forward current divided by reverse current and vice versa) which oscillates with the 
same period as the current density and which is considerably enhanced at its maxima. This 
oscillation of the asymmetry ratio as a function of b is shown in Fig. 2(b). Hence, choosing 
the appropriate parameters of the layer thicknesses and the barrier heights, the asymmetry 
can be enhanced significantly (more than one order of magnitude) , leading to the I-V char- 
acteristics presented in Fig. 1(b) reminiscent of a diode. Similar characteristics were found 
for double barrier structures with the same barrier heights but different barrier thickness. 
In contrast, for a symmetric double barrier structure (LV 2 = U4 and a = c), the phase shift 
in the forward and reverse current densities is zero, resulting in a symmetric I-V curve for 
positive and negative applied voltage. Consequently the diode behavior is lost. 

Replacing in the asymmetric double barrier structure the outer layers Mi and M5 by 
ferromagnetic metals, it is found that the diode efficiency can be controlled by an applied 
magnetic field. Furthermore, it is found that the tunneling magnetoresistance (TMR) ratio 
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itself depends strongly on the direction of the current yielding a high asymmetry ratio. The 
current density of such a magnetic double barrier structure was calculated for the case of 
identical magnetic layers using Fermi wave vectors k E = 1.09 A" 1 and kg = 0.42 A -1 which 



correspond to the spin-split free-electron-like ci-electron bands of Fe [23]. The parameters 
for the other layers are the same as for the non-magnetic case discussed above. The TMR 
ratio is defined as 

TMR 



Jap) 



y ju 

where jpt ap \ is the current with spin a for parallel (antiparallel) alignment of the magneti- 
zations in the magnetic layers. In Fig. 3(a) this TMR ratio is presented for an asymmetric 
structure with different barrier heights as a function of the thickness b. The TMR ratio 
reflects the oscillations of the current density as well as the phase shift between the forward 
and reverse bias (compare Fig. 2). It follows that the TMR ratio is also very asymmetric and 
the corresponding asymmetry ratio can reach values up to 200 (Fig. 3(b)) at the appropriate 
thickness b. 

This TMR asymmetry leads to a "magnetically controlled" diode, whose blocking effi- 
ciency can be varied by a magnetic field. To illustrate this, we calculate the value of relative 
magnetoasymmetry (RMA) 

( •dir(inv) I ■inv(dir)\ / jdir(inv) J •inv{dir)\ 

RMA = ^ /Jp -p T . \ U Z- s P (3) 

/ -dirxinv) I ■inv(dir)\ v ' 

\Jo-P I Jap ) 

where j p hr ( mv ) and j a h p r<yinv " > are the total forward (reverse) current for parallel and antiparallel 
configuration of the magnetization, respectively. In Fig. 4 the dependence of the RMA is 
shown as a function of b for the double barrier structure of Fig. 3. At the maxima the 
magnitude of the asymmetry ratio can be doubled by applying a magnetic field. 

It is noted, that much higher asymmetries can be obtained when the scattering in the 
metallic layers is weak. However, in this case the resonance peaks become narrow and 
are therefore more difficult to detect experimentally. More realistic is the case of strong 
scattering which leads to a broadening of the resonance peaks and consequently reduces 
the asymmetry ratios, with the advantage of being easier to detect experimentally. More 
critical for possible experimental observation of the predicted diode behaviour are barrier and 
metallic spacer thickness fluctuations. Further calculations show that thickness fluctuations 
which do not extend over two atomic layers preserve the quantum well states. Details of 
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FIG. 4: The relative magnetoasymmetry (RMA) as a function of the thickness b for the same 
asymmetric double barrier structure as for Fig. 3. The RMA is calculated from expression (|3|) with 
the current asymmetry ratio defined as {jf^l ( sond line) and (j^a P )/jp% P )) (dashed line). 

this work will be published elsewhere. As a final point we would like to underline that the 
general expressions (0) describes properly both ballistic and diffuse tunneling regime through 
the system. The transition from one regime to another can be retrieved if the scattering 
in the middle metallic layer M3 is zero (I3 — > 00). In this case, the problem can be solved 
exactly and the analytic expressions in (|2|) yield that all tunneling probabilities are equal 
(Dfj = D u for i,j = 2,4) so that the current density can be written as 

f = 32 = 3l = ^ / dE[f(E) - f(E + eV ext )} J D°{E, k)k<Ik 

and the condition of constant current across the whole double barrier structure is automati- 
cally fulfilled. This result means that in the case of an ideal structure without scattering, the 
purely quantum-mechanical problem of electron tunneling through the entire double barrier 
structure is solved exactly describing the direct coherent process. In this case, the voltage 
drop in each barrier is proportional to its thickness as in refs. Q and 0. This remarkable 
result is important since it directly shows how the scattering in the middle layer destroys 
the direct ballistic process so that the tunneling across the structure is not described any- 
more by a single transmission matrix but by resistors in series. Moreover, it means that the 



S 



expression (Q) can be written in the form: 

J dE[f(E)-f(E + eV ext )} J D% m) (E,K)KdK 

(4) 



•?2(4) 



+^ | d£[f - f(E 2 )} J [D° m) (E, k) - D° 2m) {E, k) 



where E\ = E, E 2 = E + eV 2 for j 2 and E\ = E + eV 2 , E 2 = E + eV ext for j 4 . The first term 
in ([|) describes the direct ballistic process across the structure and the second one describes 
local processes which appear only in presence of the scattering in the middle metallic layer. 

In conclusion, we presented a quantum theory of the tunnel magnetoresistance in mag- 
netic double barrier structures and predicted a strong "diode"-like behaviour for the I-V 
characteristics and the TMR in the case of asymmetric barriers. It was shown that the 
asymmetry ratio can be controlled by an applied magnetic field. This phenomenon is due 
to the different phase shift of the quantum well states in the middle metal layer under for- 
ward and reverse applied voltage. This structure should have an important application as a 
blocking device in a MRAM. 
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